Core-Level X-Ray Photoemission Satellites in Ruthenates: 
A New Mechanism Revealing the Mott Transition 
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Ru 3d core-level x-ray photoemission spectra of various ruthenates are examined. They show 
in general two-peak structures, which can be assigned as the screened and unscreened peaks. The 
screened peak is absent in a Mott insulator, but develops into a main peak in the metallic regime. 
This spectral behavior is well explained by the dynamical mean-field theory calculation for the single- 
band Hubbard model with on-site core-hole potential using the exact diagonalization method. The 
new mechanism of the core-level photoemission satellite can be utilized to reveal the Mott transition 
phenomenon in various strongly correlated electron systems, especially in nano-scale devices and 
phase-separated materials. 

PACS numbers: 71.10.Fd,71.30.+h,79.60.-i 



Core-level x-ray photoemission spectroscopy (XPS) 
has long been a powerful tool to investigate the chemical 
environment of solids and molecules . Usually the vari- 
ation of the binding energy depending on the chemical en- 
vironment ( "chemical shift" ) is utilized for this purpose, 
but satellite structures arising from the Coulomb attrac- 
tion between valence electrons and a core hole created 
in photoemission can be a fruitful source of information, 
especially in strongly correlated systems such as rare- 
earth and transition-metal compounds (TMCs) 0, 3. 
A good example is the high- Tc and related cuprates ^, 
where the satellite structures of Cu 2p core-level XPS 
have been instrumental to confirm the charge-transfer 
nature of the insulating gap in the parent compounds 
which were traditionally classified as Mott insulators , 
and determine basic parameters such as the d-d Coulomb 
' correlation energy and the charge-transfer energy. 

Such a strong electron correlation effect has been con- 
.sidered to be weak in 4d TMCs because Ad orbitals are 
fairly delocalized, and RuO?, for example, is traditionally 
.classified as a band metal p- However, recent studies re- 
•vealed that many ruthenates such as Ca2-2:Sr2:Ru04 
' and pyrochlores |^ have various interesting properties 
■related with the correlation effect among Ru Ad electrons. 
I Ru 3d core-level XPS spectra also show some hint of the 
strongly correlated nature of Ru Ad electrons. As shown 
in Fig. 1, the Ru 3d spin-orbit doublet, of which split- 
ting is about 4 eV, consists of roughly two components, 
"screened" (denoted as "s") and "unscreene d" ( denoted 
as "u" ) peaks in various ruthenates 0, includ- 
ing Ru02 where their origin has been the subject of a 
long debate ES [H E [H [11 . The peak "s", located 
at about 2 eV lower binding energy than that of the 
rather broad unscreened peak "u" , is absent in an in- 
sulator Y2RU2O7 and grows as the system becomes more 
metallic. Hence Cox et al. suggested that the peak "s" 
might originate from the screening of a core hole by a 
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FIG. 1: Ru 3d XPS spectra of RUO2 (after Bi2Ru207 
(after JJ), Sr2Ru04 (after O), SrRuOg, CaRuOg, 
Y2RU2O7 (after 11]). "s" and "u" denote screened and un- 
screened peaks, respectively, following TT]. All the spectra 
are ahgned by the unscreened-peak positions of Ru 3d5/2. 



quasiparticle at On the other hand, we may 

attribute it to the charge-transfer mechanism as in 3d 
TMCs, but the energy separation between the main and 
satellite peaks is much smaller than that of a Pd 3d spec- 
trum in PdO [isj . while we expect larger separation ac- 
cording to the chemical trend jSj . Quite recently, Okada 
141 tried to explain the two-peak structure in Sr2Ru04 
1^ by the "non-local screening" effect 0| in addition to 
the charge-transfer mechanism. He could not, however, 
fit the strong intensity variation and the little change of 
the energy separation in ruthenate series with reasonable 
values of parameters. 

In this Letter, we propose a new mechanism for the 
two-peak structures in the Ru 3d core-level XPS spectra 
of ruthenates based on a single-band Hubbard model. 
The core-level spectra of a single-band Hubbard model 
with core-hole potential are calculated by the dynami- 
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cal mean-field theory (DMFT) 17]. The DMFT can de- 
scribe successfully the Mott transition with the change 
of ratio of the Coulomb interacion U between conduc- 
tion electrons to the band width W, and it should be 
suitable for a local problem like core-level XPS despite 
its shortcomings due to "local" electron self-energy |l8|. 
In this picture, the "screened" peak originates from the 
screening of a core-hole by quasiparticles on the Fermi 
surface. We found that the calculation can reproduce 
the spectral behavior shown in Fig. 1 just by changing 
the correlation strength with reasonable parameter val- 
ues. We also made a systematic study of the Ru 3d 
XPS spectra of a bandwidth-controlled Mott-Hubbard 
system Y2_2:Bia;Ru207 ;9| as a function of composition 
to strengthen our conclusion. 

We start from the single-band Hubbard model 



(1) 
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where the notation is standard, which is exactly mapped 
in infinite dimensions on to the sing le- impurity Anderson 
model according to the DMFT [l7| 



1=2,17 a 

■J2Viicld, + n.c.). (2) 
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In order to describe a core-level spectrum, we should add 
a term due to a core hole for XPS final states after getting 
the converged model parameters. 



^coro = (e/i - Q^dld^)h^h, 



(3) 



where the operator {h) creates (destroys) a core hole 
on the "impurity" site, eh is a core-hole energy, and Q (= 
1.25U, for simplicity in all calculations) is the Coulomb 
interaction between the core-hole and "impurity-site" 
conduction electrons. In order to solve the model with 
the semi-elliptical shape for the conduction-band density 
of states, we followed the exact diagonalization method 
proposed by CafFarel and Krauth [Igj . 

Figure 2 shows core-level and valence-band spectra of 
the ten-site {Ns = 10) half-filled model obtained by the 
modified Lanczos method changing the ratio W/U from a 
Mott insulator regime to a good metal regime to show the 
behavior of spectral weights in the bandwidth-controlled 
metal-insulator transition. For this particle-hole sym- 
metric model, ed is chosen to be —U/2, then the chem- 
ical potential fj, is identically zero. The core-hole en- 
ergy Ch is set to Q. The core- level spectra are mainly 
composed of two peaks (minor peaks are due to the fi- 
nite size) and their behavior with W/U is quite similar 
to that of the valence-band spectra and, moreover, the 
intensity of the "screened peak" is found to be nearly 
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FIG. 2: (a) Core-level and (b) valence-band spectra of the 
half-filled Hubbard model with Ns = 10 calculated by the 
DMFT varying W/U. Solid lines are guidelines obtained by 
broadening dotted lines with a Gaussian of 0.25 full width at 
half maximum to remove discreteness due to the finite size. 



proportional to the quasiparticle weight Z calculated by 

(1 — ^^^^^r^) "'^ L^o' ^-'^i'^-'^ implies that the two-peak 
structure is strongly related with the Mott transition. In 
a Mott insulator {W/U = 0.65) we can see a broad single 
peak, and as the coherent peak grows with the increase 
of W in the valence-band spectra, a sharp shake-down 
satellite structure appears in the core-level spectra in the 
positive energy side and eventually becomes a main peak 
with the asymmetric low-energy tail. Around W/U = 2.0 
(top-most figure), both the core-level and valence-band 
spectra remind of the famous 6 eV satellite in Ni metal, 
which is explained by the Mahan-Nozieres-DeDominicis 
model with strong core-hole potential for the former and 
by the Hubbard model for the latter Now, both 

spectra are explained in a single model thanks to the 
DMFT based on the exact diagonalization, and it should 
be applicable to other strongly correlated systems show- 
ing the Mott transition. 

In order to investigate in detail the origin of the two- 
peak structure in the core-level spectra, we obtained all 
the eigenvalues and eigenvectors of the initial and final 
states of a four-site model by explicit exact diagonaliza- 
tion. Figure 3 shows the calculated core-level spectra for 
the case of (a) W/U = 0.75 and (b) 1.5. Here we see again 
that for a narrow band case (a) only the unscreened peak 
is dominant, while for the wide band (b) both screened 
and unscreened peaks appear in the core-level spectra. 
To understand the valence-electron configuration of these 
different core-level final states, we plot at the bottom of 
the figure the bar diagram showing the weights of various 
valence-electron configurations contributing to each core- 
level peak. For the four-site model, we should consider 
three energy levels in the conduction band representing 
a lower Hubbard band (LHB), a coherent peak, and an 
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upper Hubbard band, along with the "impurity" level at 
the core-hole site in the presence of a core-hole. Then 
there are three major valence-electron configurations in 
the XPS final state, i.e. (= d^L^C^ for iV^ = 4), d^C, 
and d^L, where d denotes the "impurity" level and C (L) 
denote the coherent peak (LHB). These valence-electron 
configurations arc shown schematically in Fig. 4. It can 
be easily noticed from the bar diagram at the bottom 
of Fig. 3 that the unscreened peak originates from the 
configuration d^, while the screened peak from the con- 
figurations d^C and d^L. 

With these configurations, the behavior of two-peak 
structure with the change of W/U is simply understood. 
When W/U is so small that the system becomes a Mott 
insulator, the configuration d^C is absent both in the 
ground and the final states, so the unscreened peak, of 
which main configuration is d^, becomes dominant, and 
the screened peak by the configuration d^L is located 
near the unscreened peak (the separation is Q — U in the 
limit W/U 0). Since the configuration d^L represents a 
single-hole state in the LHB, the core- level peak of a Mott 
insulator has extra lincwidth broadening reflecting the 
LHB shape. When W/U is close to the critical value and 
the system becomes a metal, the configuration d^C is now 
available to take part in the screening of the core hole, 
which makes a sharp shake-down satellite peak, whose 
energy is about Q — U/2 lower than that of the unscreened 
peak and its intensity increases with the weight of the 
coherent peak. 

To test the validity of the present model, we 
made a systematic study of the Ru 3d core-levels in 
Y2_a:Bi3:Ru207 systcm. This system is known to be a 
bandwidth-controlled Mott-Hubbard system from the- 
oretical calculation j^, transport properties ^ and 
valence-band photoelectron spectra , which shows the 
metal-insulator transition around a; = 0.9 0- We took 
Ru 3d XPS spectra of polycrystalline Y2-a;Bia;Ru207 




FIG. 4: Schematic diagrams of three major valence-electron 
configurations contributing to the core-level spectra. Note 
that the "impurity" level at the core-hole site is lowered by 
the attractive Coulomb potential Q. 

samples (x = 0.0, 0.4, 1.0, 1.6, and 2.0) using Mg Ka 
source after annealing four hours under the oxygen par- 
tial pressure of 10"^ Torr at 900 K. The anneahng pro- 
cess, which removes almost all of the carbon contami- 
nations in the sample, is essential because the binding 
energy of the Ru 3d core-level is close to that of C Is 
peak. 

Ru 3c? XPS spectra (solid lines) of Y2_a;Bia;Ru207 thus 
obtained are shown in Fig. 5 (a) after removing the inelas- 
tic background and the small residual C Is peak contribu- 
tion by deconvolution using Doniach-Sunsic line shapes. 
We can see the systematic behavior that the weight of 
the screened peak becomes larger as the Bi concentra- 
tion increases. This tendency is strongly correlated with 
the transport properties Q and valence-band photoelec- 
tron spectra [2l| . Using the present model, the core-level 
spectra of Y2-2:Bi2:Ru207 series could be reproduced by 
changing only W with fixed U value (=1.7 eV). We also 
fit the spectra of other ruthenates shown in Fig. 1 by 
freeing U, and Table 1 summarizes all the parameter val- 
ues. After proper broadening of calculation results by 
a Voigt function to simulate the core-hole lifetime and 
the experimental resolution, final results (dotted lines) 
are plotted over the experimental spectra in Fig. 5 (a) 
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FIG. 3: Core-level spectra in the four-site model for (a) 
W/U — 0.75 and (b) 1.5. The calculated spectra are broad- 
ened by a Lorentzian of 0.25 full width at half maximum. 
The bar diagram at the bottom indicates the weights of dif- 
ferent valence-electron configurations in each XPS final state. 
The three valence-electron configurations d^, d^C, and d^L 
are shown schematically in Fig 3. 
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FIG. 5: Comparison of Ru 3d XPS spectra (solid lines) of 
(a) Y2-2!Bia;Ru207 and (b) other ruthenates with the present 
model calculations (dotted lines). See text for details. 
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Y2_a:Bia;Ru2 07 and (b) other ruthenates. The fitting re- 
sults are quite satisfactory and all the spectral behaviors 
are well reproduced. Rather surprisingly, Ru02, which 
was believed to be well described by band calculations 
01 , should be reconsidered as a strongly correlated metal 
according to our fitting result {W/U = 2.0). 

Although the present model does not include such 
terms as 4(i-orbital degeneracy and also the Ru ^dt^^ 
band shape is far from semi-elliptical [j^ E^j ). it is 
quite successful in describing the two-peak structure of 
the core-level spectra, which confirms that the Mott- 
Hubbard mechanism is the most important factor. Nev- 
ertheless, the role of O 2p bands could not be ruled 
out in two respects: (1) The change of VF- value for 
Y2_2:Bi:cRu207 is rather larger than expected 0, E^ . 
while W values for other ruthenates are comparable to 
those of band calculations 23] . Since we are dealing with 
an effective Hubbard model, in which O 2p bands are in- 
tegrated out, the U value should change in real materi- 
als depending on the Ru Ad-0 2p hybridization strength, 
which would result in the reduction of the W variation. 
(2) We may expect that a similar two-peak structure 
should appear in other strongly correlated TMCs. Actu- 
ally, there is some hint of it in Ti203 Q and V2O3 
But the screened-peak intensity is rather small in these 
systems compared with the coherent peak of the valence 
band |2^ . which may be explained by the competition 
between the charge-transfer screening by the O 2p elec- 
tron and the screening by quasiparticles at Ep. Since 
the quasiparticle screening takes place through mediat- 
ing O 2p orbitals, it is probably more rapidly suppressed 
than the charge-transfer screening as the hybridization 
strength becomes weak. Further experimental and theo- 
retical studies are necessary to clarify this issue. 

In summary, we have elucidated the origin of two-peak 
structures in the Ru 3d core-level XPS spectra of various 
ruthenates as due to the different screening mechanisms 
in the Mott-Hubbard picture of the metal-insulator tran- 
sition. This mechanism should be applicable to other 
TMCs as well, hence core-level XPS can be utilized to 
distinguish metallic from insulating regions in phase- 
separated materials j2^ and TMC-based nano-clcctronic 
devices when used in conjuction with spectro- microscopic 
tools such as scanning photoemission microscopy. 

This work was supported in part by the Korea Sci- 



TABLE I: Parameter values (in eV) obtained from fitting Ru 
3d XPS spectra of ruthenates. 
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